Electrically controlled superconducting states at the heterointerface SrTiOa/LaAlOs 
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We study the symmetry of Cooper pair in a two-dimensional Hubbard model with the Rashba- 
type spin-orbit interaction as a minimal model of electron gas generated at a heterointerface of 
SrTiOa/LaAlOs. Solving the Eliashberg equation based on the third-order perturbation theory, we 
find that the gap function consists of the mixing of the spin-singlet dxy-wsLve component and the 
spin-triplet {px ±ipy)-wave one due to the broken inversion symmetry originating from the Rashba- 
type spin-orbit interaction. The ratio of the d-wave and the p-wave component continuously changes 
with the carrier concentration. We propose that the pairing symmetry is controlled by tuning the 
gate voltage. 
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Recent development of technology of epitaxial growth 
makes it possible to fabricate heterointerface between two 
different transition-metal oxidesJ^ The discovery of the 
generation of two-dimensional electron gas (2DEG) at 
the 71- type heterointerface SrTiOa/LaAlOa attracts much 
intereslj^ since both SrTiOa and LaAlOa are band insu- 
lators in the bulk material. It is also noted that this 
2DEG shows superconductivity with transition tempera- 
ture Tc = 200 mK.— This superconductivity differs from 
the superconductivity of bulk SrTiOa with Tc = 300 mK 
in the carrier concentration and the dimensionality of 
mobile electronsi^ The discovery of superconductivity is 
significant since it is promising to control the supercon- 
ducting states electrically by the applied voltage. Con- 
trolling the number of carrier by applied voltage in super- 
conductors is one of the intriguing issues to understand 
the mechanism of the superconductivity. In cuprates. a 
superconducting region is located next to an antiferro- 
magnetic phase in an n-T phase diagram. However, in 
cuprate superconductors, it is not easy to tune the num- 
ber of carrier by external fields such as the magnetic field 
or the gate voltage since the number of carrier depends 
on material compositions. On the other hand, the carrier 
number of 2DEG at the heterointerface SrTiOs/LaAlOs 
can be tuned by the gate voltage since the carrier is gen- 
erated by the polar and unpolar nature of LaAlOa and 
SrTiOa, respectively Actually, electrostatically tuned 
superconductor-metal-insulator transition was found to 
be accessible ji^i^ 

When the electron correlations (Coulomb interaction) 
are responsible to the superconductivity, the pairing sym- 
metry depends on the number of carrier: The spin-triplet 
p-wave pairing is favored at low density )^°i^^'^^i^^'^^ while 
on the other hand, the spin-singlet d2;2_j^2-wave pairing 
is favored near half filling.— Thus, at the heterointer- 
face of oxides with large Coulomb interaction, the phase 
transition between the spin-singlet and spin-triplet states 
can be observed by the applied voltage. 

Another characteristic feature at the heterointerface 
is lack of the inversion symmetry along the direction 
perpendicular to the interface (z-axis), which induces 
the asymmetric spin-orbit interaction called the Rashba- 
typc spin-orbit interaction (RSOI) This feature is com- 



mon to non-centrosymmetric superconductors such as 
CePtgSi, CcRhSia, and CelrSig ^^i^i^o In these super- 
conductors, the admixture of even- and odd-parity pair- 
ings has been suggested theoretically.— i^iS^ii^ Thus, we 
expect the mixing of parity also occurs in the present 
system. In this case, the change in the pairing symme- 
try under applied voltage is not a phase transition but a 
crossover between different parities. Nevertheless, volt- 
age control of the pairing symmetry in 2DEG generated 
at the oxides interface could be a useful tool to study the 
carrier concentration dependence of the superconducting 
state. Thus, it is very timely to reveal the pairing sym- 
metry of Cooper pair in 2DEG generated at the oxide 
interface. 

In this Rapid Communication, we study the pair- 
ing state of superconductivity at the heterointerface of 
SrTiOa/LaAlOa by using the Hubbard model with the 
RSOI. Since both SrTiOa and LaAlOs are band insula- 
tors in the bulk, the number of carrier in 2DEG is far 
from the half filling. In this case, the inherent strongly 
correlated nature of the titanium oxides may not emerge. 
Therefore, we solve the Eliashberg equation based on per- 
turbation theory up to the third-order expansion of the 
Coulomb potential. We find that the pairing state in this 
system is the admixture of spin-singlet even-parity {dxy- 
wave) pairing and spin-triplet odd-parity {px ± ipy-wave) 
pairing, which arises from broken inversion symmetry at 
the heterointerface. The ratio of these two components 
continuously changes with the number of carrier. 

The mobile electrons at the heterointerface are mainly 
introduced to 3d orbitals of Ti'^+ ions. Considering 
the crystalline electric field by an ionic model, the two- 
dimensional dxy orbital has the lowest orbital energy level 
in 3d orbitalsi^ In fact, it has been confirmed that the 
dxy electrons compose the 2DEG at the interface by x- 
ray absorption spectroscopy with the linearly polarized 
light 1^ Several results in density- functional calculations 
agree with this experimental onci^i^ Thus, we consider 
a two-dimensional Hubbard model with the RSOI as a 
minimal model of 2DEG generated at the heterointerface 
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of SrTiOa/LaAlOs. The Hamiltonian is given as 
H = -t (cLcjT + /I.C.) +^C/n,|n,x 



A 
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(1) 



where Ci„{cl^) is an annihilation (a creation) operator 

of an electron with spin <t at site i, and riia = cl^Cia- 
denotes the set of the nearest neighbor sites and 
t is the transfer integral. The third term is the RSOI, 
where A is the magnitude of Rashba field and & are the 
Pauli matrices. The vector g{k), which satisfies the rela- 
tion g{k) — ~g{~k), induces the breakdown of inversion 
symmetry. We adopt g{k) ~ {~Vy{k),Vx{k),0)/t, with 
Vx,y{k) = dek/dkx^y = 2tsa\{kx.y). The bare Green's 
function is given by the following 2x2 matrix in spin 
space, 



G{k, iSji) — 



^ f G'i^{k,ien) G'|x(fc,ie„) 



£k)I+^g{k) ■ (T 
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where I is the unit matrix, and e„ ~ (2n — 1)t:T is the 
Matsubara frequency for fermions. 

The effective pairing interaction Vaj^^^{k,k') in this 
study is given by the perturbation expansion up to the 
third order with respect to U, 

V,,,Ak,k') =V^l^ik,k') + V^;:'^^-{k,k'), (3) 

V^l^ik, k') =U + U\{k + fc') + U^^ik + fc-') 

+ U^xHk-k'), (4) 

V^;_T^ik, fc') =2[/3Rc^G(fc - q)x{q)G{k' - q) 
q 

-2U^R(Y2G{k + q)(p{q)G{-k' +q), (5) 



V^^^ik,k')^-U\{k-k'), 



(6) 



Vertex 



(fc, fc') =2[/3Rc^G(fc + q)xiq)G{k' + q) 



+ 2C/3Re^G(-fc + q)<j)iq)Gi-k' + q), (7) 
q 

where G(fc) = GTT(fc) = Gxi(fc), x{q) = ~T/Nj:,G{q + 
k)G{k), and 0((j) = -T/NJ2kG{q - k)G{k). fc = 
(fc, {2n + 1)ttT) and q = {q, 2nniT) are the short no- 
tation of momentum and Matsubara frequency. Here, we 
have dropped the terms including the off-diagonal part 
of G(fc) for simplicity. We divide Vcria2{kik') into two 
parts yJ^^P^(fc,fc') and V^^ZT^'ik.k') for the later discus- 
sions. V"j^^^(fc, fc') are the term included in the RPA, and 
V^^^l°^{k, fc') are the contributions from the other terms. 

In order to study the pairing instabilities, we solve the 
linearized Eliashberg equation by a power method, 

T 



N 



y ] Vai(72 {k, fc )Faia2 {k ), 



(8) 



i(T2(^) — ^ ^ Ger^(73 (fc) Aa-3g-4 (fc)G,,,,(-fc). (9) 



Solving Eqs. ^ and ([9]) self-consistently, we obtain the 
eigenvalue a of the Eliashberg equation and pairing func- 
tion Acri(j2(fc)- The superconducting transition tempera- 
ture Tc is identified as the temperature with a = I. Thus, 
the pairing symmetry with largest value of a is the most 
dominant. Here, the indices of spin of the right-hand 
side and the left-hand side in Eq. ([S]) are the same since 
we have dropped the spin-flip scatterings in the effec- 
tive interaction. Therefore, the admixture of pairings is 
not induced by the effective interaction V. However, the 
spin-flip processes in Eq. ^ induce the admixture of 
pairings with different parities. 

The spin-singlet and the spin-triplet components with 
^2 = are extracted by A||(fc) ± Aj|(fc) where -1-1 and 
— 1 correspond to triplet and singlet ones, respectively. 
The spin-triplet components with Sz — ±1 are given 
by A||(fc) and Ajj^(fc), respectively. Due to the broken 
inversion symmetry, the components of odd-frequency 
gap function, which are extracted by {A^j^^^{k,ien) — 
Acri(j2 (fc, — jeTi)}/2, become finite in addition to even- 
frequency onci^ However, we do not mention the odd- 
frequency component hereafter since the magnitude of 
the odd-frequency component is much smaller than that 
of the even-frequency one. We choose t = 1 for the unit of 
energy. We take 128 x 128 fc meshes and 2048 Matsubara 
frequencies. 

Att(fc,ie„) - Ai^Xk,i£„) A^4(fc,i£„) - Ai,^ik,ie„) 

-0.0008 




FIG. 1: (Color online) The fc-dependence of the spin-singlet 
gap function Af i{k, ieo) — A if{k, ieo) and the spin-triplet gap 
function An(fe,ieo) - Au(fe,ieo) for T/t = 0.008, X/t = 0.1, 
U/t — 3 and n — 0.3, with eo = ttT. The dotted (green) and 
solid (black) lines denote the nodal lines and Fermi surfaces, 
respectively. 



First, we look at the gap function. In the absence of 
the RSOI, the spin state of Cooper pair can be classi- 
fied into spin-singlet or spin-triplet states. For n < 0.5, 
the symmetry of the gap function with the largest val- 
ues of a is dxy-'wave singlet or p-wave triplet pairing. By 
introducing the RSOI, the pairing functions with differ- 
ent parity are mixed due to the breakdown of inversion 
symmetry, since Ga,-a{k) becomes finite in Eq. ([9]). For 
this reason, the dxy-wave singlet pairing and the p-wave 
triplet pairing with 5*2 = ±1 coexist. Indeed, as shown in 
Fig. [U we have confirmed that the pairing symmetry of 
the A|j(fc) — A|-f (fc) is dxy wave, and that of A||(fc, ie„) 
and Ax|(fc, ze„) are (px + iPyY and {—px + ipy) waves. 
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respectively. Here, the component of p-wave triplet pair- 
ing with Sz = does not exist, since the RSOI induces 
the mixing between the pairings with different parity of 
momentum and different values of 
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FIG. 3: (Color online) The ratio of the singlet and the triplet 
components for T/t — 0.008 and A/f = 0.1 as a function of 
U. 



FIG. 2: (Color online) Eigenvalues of Eliashberg equation 
against n for T/t = 0.008 and U/t = 4. 



Next, in Fig. ^ we show the n-dependence of a at 
U/t = 4 and T/t = 0.008. At first, we concentrate on 
the case without the RSOI. The value of a for the dxy- 
wave pairing is the largest at n ~ 0.5, while on the other 
hand that for the p-wave pairing becomes the largest for 
n ^ 0.4. This result can be understood as follows. Since 
V^l^ik,k') is the repulsive force {V > 0), V^l^{k,k') 
favors the singlet pairing whose gap function changes the 
sign through Q, where Q is the wave vector at which the 
magnitude of x(Q, itOm) is large. Since the magnitude of 
V^^^{k, k') increases with n, the value of a for the d^y- 
wave pairing also increases. On the other hand, the calcu- 
lated V^°^^'^^{k, k') is attractive for the p-wave pairing for 
n < 0.5. Since the relative contribution of V^^^l'^^{k, k') 
in Va^a^i^^k') increases with decreasing n, the p-wave 
pairing becomes dominant at low carrier concentration. 
This is consistent with the previous studies using pcrtur- 
bative approaches including vertex terms.— liii^ii^ No- 
tice that the values of a for d2;2_j,2-wave pairing are 
smaller than those for d^-y-wave pairing for n ^ 0.5 since 
V^'f1j°^{k, k') suppresses it. In the present calculation, 
the dx^-y^-wave pairing becomes dominant for n ^ 0.6 
as noted in the previous studies in the context of high-Tc 
cuprateSfi^ii^ 

Next, we focus on the case oi X/t ~ 0.1. In the presence 
of the RSOI, breakdown of inversion symmetry induces 
the mixing of the even- (d-wavc) and the odd-parity {p- 
wave) pairing. For the small magnitude oi n (n ^ 0.4), 
the main component of the gap function with largest a is 
the spin-triplet p-wave one. This is because Vrr^a^i^ik') 
is attractive for p-wave pairing. On the other hand, for 
n ^ 0.4, V becomes attractive for dxy-wave pairing and 
therefore, the magnitude of the d-wave component be- 
comes dominant. Thus, the behavior of n-depcndence 
of a is similar to that for the d^^y-wave and the p-wavc 
pairings without the RSOI for n < 0.4 and n > 0.4, re- 




FIG. 4: (Color online) n-U phase diagram of the symmetry 
of gap function at T/t = 0.008. The (red) dashed and (red) 
dotted lines denote the lines where the corresponding values 
of K equal to and ±0.4, respectively. The (black) solid line 
denotes the boundary between spin-singlet and spin-triplet 
phase in the absence of the RSOI. 



We have confirmed that the magnitude of the spin- 
singlet (da;j,-wave) and spin-triplet (p-wave) components 
changes with n. Next, in Fig. [31 the ratio of singlet 
and triplet components k = (|At| — |As|)/(|At| -I- \As\) is 
plotted as a function of U, where |As_t| are defined as the 
absolute value of the maximum of the singlet and triplet 
components, respectively. The values of k becomes unity 
for purely triplet case and minus unity for purely singlet 
case. In the presence of the RSOI, k increases with U 



since the magnitude of V^^'^^^°^{k, k') dominates over that 
of V^^^[k^k'). This is because the third-order term in 
7^aT''{k.k') is larger than that in V^^J^{k,k') for n < 



0.5. In the absence of the RSOI, the values of k change 
from —1 (pure singlet) to 1 (pure triplet) abruptly. On 
the other hand, the value of k varies continuously as a 
function of n in the presence of the RSOI. 

Here, we discuss the possibility of the crossover be- 
tween the spin-singlet and the spin-triplet states at the 
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actual heterointcrfacc SrTiOs/LaAlOs. In Fig. [4l we 
show the n-U phase diagram of the pairing states. The 
dashed hne denotes the hne where the magnitude of the 
spin-singlet and the spin-triplet components are the same 
(k = 0). This line corresponds to the dashed line in Fig. 
[3l In both Figs. [3] and [H as a guide of eyes, we also 
draw the dotted lines with k ~ ±0.4, where the corre- 
sponding values of the ratio |At| : |As| are 7:3 and 3:7 
for plus and minus sign, respectively. We also depict the 
boundary between spin-singlet and spin-triplet phase in 
the absence of the RSOI by the solid line. In 3c? orbitals, 
the value of U/t is popularly thought to be 3-4, which 
is approximately half of the band width. In these values 
of U/t, the sign of k changes at around n = 0.3-0.5. If 
the number of carrier is controlled around these values 
by the applied gate voltage, the crossover behavior be- 
tween the spin-singlct and the spin-triplet states can be 
observed. Note that n = 0.5 is an ideal value for intrin- 
sic doping, which originates from the polar and unpolar 
nature of LaAlOa and SrTiOa, respectively<^ In order to 
observe this crossover, we need precise control of the car- 
rier number and/or the large magnitude of the RSOI. 
If the magnitude of the RSOI is small, the value of k 
varies abruptly with n like the case without the RSOI. 
Therefore, the crossover behavior in a small magnitude 
of the RSOI is similar to a phase transition between the 



spin-singlet and spin-triplet states. Since the qualitative 
results in the present calculation do not change with the 
values of U, A and T, the crossover between the spin- 
singlet and spin-triplet states might be observed at the 
heterointerface SrTiOs/LaAlOa or the related materials. 
The direct determination of the superconducting energy 
gap from bulk property is not easy since the transition 
temperature is very small in the present system. One of 
the possible ways is to detect spin current via Andreev 
bound state. We can expect the enhancement of the re- 
sulting spin current when the magnitude of the triplet 
component of the pair potential is larger than that of the 
singlet onci^i^ 

In this Rapid Communication, we have studied 
the pairing symmetry in two-dimensional Hubbard 
model with the RSOI considering the heterointerface 
of SrTiOs/LaAlOa. Solving the Eliashberg equation 
based on the third-order perturbation theory, we have 
found that the gap function consists of the mixing of 
the spin-singlet d^jy-wave component and the spin-triplet 
(px ± zpy)-wave one because of the broken inversion sym- 
metry originating from the RSOI. The ratio of the d- 
wavc and the p-wave components continuously changes 
with the carrier concentration through the change in the 
effective pairing interaction. 
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